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f3 ■ Abstract. 

We report on investigations of the chiral and deconfinement aspects of the finite temperature 
*-M transition in 2+1 flavor QCD using the Highly Improved Staggered Quark (HISQ) action on 

lattices with temporal extent Nt = 6 and Nt = 8. We have performed the calculations for 
physical values of the strange quark mass rris and the light quark masses rrii — 0.2ms and 
0.05ms. Several finite temperature observables, including the renormalized Polyakov loop, the 
renormalized chiral condensate and the chiral susceptibility have been calculated. We also study 
^ , the fluctuations and correlations of different conserved charges as well as the trace anomaly 

at flnite temperature. We compare our flndings with previous calculations that use different 
improved staggered fermion formulations: asqtad, p4 and stout. 
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^ ■ 1. Introduction 

Quantum Chromodynamics (QCD) is the theory of strong interactions that successfully describes 
a wealth of experimentally observed phenomena. At high energies the coupling of the theory 
^ ■ runs to zero [l], making interactions between the constituents of the theory (quarks and gluons) 

weaker, eventually approaching the limit of non-interacting ideal gas. 

In the region of low energies (up to a few hundred MeV) the coupling of the theory is of 
0(1) (thus justifying the name strong interactions) and the perturbation theory is not reliable. 
In this phase of the theory the degrees of freedom are hadrons, while quarks and gluons are 
confined into objects with zero net color charge. 

The outlined picture suggests that at some temperature a transition from the confined 
hadronic phase to the deconfined quark-gluon plasma (QGP) phase takes place [2j. Creating 
conditions for achieving QGP has been the subject of recent experiments at the Relativistic 
Heavy-Ion Collider (RHIC) [3j at the Brookhaven National Laboratory and is going to be the 
goal of the future heavy-ion program at the Large Hadron Collider (LHC) [4J at CERN. 

Since the physics of the hadronic phase and the transition to QGP is essentially non- 
perturbative, lattice QCD techniques have been extensively applied to its studies [5]. Putting 
fermions on the lattice leads to the infamous fermion doubling problem, and several lattice 
fermion formulations deal with it in different ways. Improved staggered fermion formulations 
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are widely used to study QCD at non-zero temperatures and densities, see e.g. Ref. [6l [7] 
for recent reviews, for, at least, two reasons: they preserve a part of the chiral symmetry of 
the continuum QCD which allows one to study the chiral aspects of the finite temperature 
transition, and are relatively inexpensive to simulate numerically because due to absence of an 
additive mass renormalization the Dirac operator is bounded from below. However, there are at 
least two problems with staggered fermion formulation. The first one is the validity of the rooting 
procedure, i.e. the way to reduce the number of tastes (unwanted unphysical degrees of freedom 
that result from fermion doubling) from four to one, and the other is breaking of the taste 
symmetry at finite lattice spacing. The discussion of the validity of rooted staggered fermions is 
presented in Refs. [8l[9]. To reduce the taste violations smeared links, i.e. weighted averages of 
different paths on the lattice that connect neighboring points, are used in the staggered Dirac 
operator and several improved staggered formulations, like p4, asqtad, stout and HISQ differ 
in the choice of the smeared gauge links. The ones in the p4 and asqtad actions are linear 
combinations of single links and different staples [10^ [TT] and therefore are not elements of 
the SU(3) group. It is known that projecting the smeared gauge fields onto the SU(3) group 
greatly improves the taste symmetry [12j. The stout action [13] and the HISQ action implement 
the projection of the smeared gauge field onto SU(3) (or simply U(3)) group and thus achieve 
better taste symmetry at a given lattice spacing. For studying QCD at high temperature it 
is important to use discretization schemes which improve the quark dispersion relation, thus 
eliminating the tree level O(a^) lattice artifacts in thermodynamic quantities. The p4 and 
asqtad actions implement this improvement by introducing 3-link terms in the staggered Dirac 
operator. In this paper we report on exploratory studies of QCD thermodynamics with the 
HISQ action which combines the removal of tree level O(a^) lattice artifacts with the addition 
of projected smeared links that greatly improve the taste symmetry. We also compare our results 
with the previous ones obtained with the asqtad, p4 and stout actions [HI [El dH dTJ UHl |T9] 
at comparable quark masses and lattice spacings as well as with the asqtad action at smaller 
quark masses and larger Nr |20j . 



2. Action and run parameters 

The Highly Improved Staggered Quark (HISQ) action developed by the HPQCD/UKQCD 
collaboration [21j reduces taste symmetry breaking and decreases the splitting between different 
pion tastes by a factor of about three compared to the asqtad action. The net result, as recent 
scaling studies show [22^ i23j. is that a HISQ ensemble at lattice spacing a has scaling violations 
comparable to ones in an asqtad ensemble at lattice spacing 2/3a. 

In this exploratory study we used the HISQ action in the fermion sector and the tree- level 
Symanzik improved gauge action without the tadpole improvement. The strange quark mass 

rus was set to its physical value adjusting the quantity J2m\ — ml = rrij^^^ ~ \/2Brn^ to the 
physical value 686.57 MeV. Two sets of ensembles have been generated along the two lines of 
constant physics (LCP): rrii = 0.2ms and mi = 0.05ms, where mi is the u and d quark mass. In 
the first set runs were performed on 16^ x 32 lattices at zero temperature and 16^ x 6 at finite 
temperature, and in the second set on 32^* and 32^ x 8 lattices, correspondingly. To estimate the 
cutoff effects we have also generated a set of high-temperature ensembles on 24^ x 6 lattices along 
the m.1 = 0.05m,s LCP. The parameters and statistics of the runs are summarized in Tabled) The 
molecular dynamics (MD) trajectories have length of 1 time unit (TU) and the measurements 
were performed every 5 TUs at zero and 10 TUs at finite temperature. Typically, at least 300 
TUs were discarded for equilibration at the beginning of the simulations. 

The lattice spacing has been determined by measuring the static quark anti-quark potential. 
As in previous studies by the MILC collaboration the static potential was calculated fixing the 
Coulomb gauge and considering temporal Wilson lines of different extent. Forming the ratio of 
these correlators and fitting them to constant plus linear Ansatz we extracted the static potential 
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Figure 1. The static potential calculated for mi = 0.2ms (left) and m,i = 0.05ms (right) in 
units of rQ. In the left figure we compare the HISQ result with the p4 result obtained at similar 
value of the lattice spacing. The dashed line on the right figure is the string potential (see text). 



V{r). We have calculated the Sommer scale vq and the related ri scale defined as 
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The static potential calculated for 0.2m.s and O-OSm^ is shown in Fig. [TJ The potential has been 
normalized to the string potential 
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string 
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+ ar, 
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at r = 1.5ro or equivalently to the value 0.91/ro at r = tq. The additive constant determined 
by this normalization is used to calculate the renormalization constant for the Polyakov loop as 
will be discussed later. There is no large cutoff dependence visible in the static potential. 
Furthermore, the static potential calculated with the HISQ action is very similar to that 
calculated with the p4 action. For the ratio of rg and ri scales we get 

ro/n = 1.459(3), mi = 0.2ms, ro/n = 1.481(6), mi = 0.05m^. (3) 

To convert from lattice units to physical units we use the value r^ = 0.469 fm determined in 



The masses of several hadrons measured on zero-temperature ensembles fall into ranges 
summarized in Table [2l They give an idea about the lines of constant physics in our calculations. 
The pion and kaon masses are most sensitive to the choice of the quark mass. We see that the 
0.2ms LCP corresponds to the lightest pion mass of about 300 MeV, while the other 0.05ms 
LCP corresponds to the lightest pion mass of about 160 MeV. The latter is very close to the 
physical value of about 140 MeV. For this reason we refer to m/ = 0.05m5 as the physical quark 
mass. From the table we see that the systematic error in the choice of the quark masses is about 
4% for the 0.2ms LCP and about 3% for the 0.05ms LCP. 

As mentioned above the lattice artifacts for the HISQ action are significantly reduced 
compared to the asqtad action. The taste violations are strongest in the pseudo-scalar meson 
sector. Therefore we studied the splitting of pseudo-scalar meson masses in the 8 different 
multiplets for m,i = 0.2ms. The quadratic splittings of the pseudo-scalar meson masses are 
independent of the quark mass to a very good approximation and therefore it is sufficient to 



Table 1. The parameters of the numerical simulations: gauge coupling, strange quark mass 
and the number of time units (TU), i.e. the number of MD trajectories for each run. Here TU, 
stands for the number of time units for zero-temperature runs, while TU, T is the number 
of time units for finite-temperature runs. The last column gives the number of time units for 
24^ X 6 runs. 
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Table 2. Ranges of masses (in MeV) of several hadrons for the two sets of ensembles. 
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study them for larger value of m,/ . The results are shown in Fig. [2] and compared to the stout 
results. In the figure we also show the masses of all the pseudo-scalar mesons as function of 
the lattice spacing for asqtad and stout actions assuming that the lightest pseudo-scalar mass 
is fixed to its physical value. For the lattice spacing corresponding to the transition region for 
temporal extent A^,- = 8 the mass of the heaviest pion is about 400 — 600 MeV for the stout and 
asqtad action even for the physical values of the light quark mass. The quadratic splittings in 
HISQ calculations are 2 to 3 times smaller than in the calculations with the asqtad action and 
also somewhat smaller than for the stout action. However, the effect of mass splitting in the 
pseudo-scalar sector is non-negligible even for HISQ. 

The smaller taste violations also result in better scaling of other hadron masses. In Fig. [3] the 
K* , p and (j) meson masses as well as the nucleon mass calculated for 0.05m, are shown. As one 
can see there is a good agreement for the vector meson masses calculated on the lattice and the 
experimental values. The nucleon mass in the HISQ calculations is about 10% larger than the 
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Figure 2. The splitting between different pion multiplets for the HISQ action at 0.2ms 
compared to the stout results shown as open symbols (left). In the right figure the non-Goldstone 
pseudo-scalar meson masses are shown as function of lattice spacings assuming that the lightest 
(Goldstone) pion mass is fixed to its physical value. The open symbols in the right figure 
correspond to pseudo-scalar mesons labeled as "1" and "3" (or, equivalently, as 7J75 and 7j7j) 
in the stout calculations. 
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Figure 3. The hadron spectrum for the 
HISQ action at rai = O.OSttt,^ compared with 
experiment. The open symbols show the 
nucleon mass for the asqtad action. For clarity 
of the plot the nucleon mass is shifted by 
150 MeV. 



experimental value, but the discrepancy between the lattice calculations and the experiment is 
largely reduced compared to the asqtad calculations. 

3. Deconfinement transition 

The deconfinement transition is usually studied using the renormalized Polyakov loop [151 (Ml EBl 
125] . It is related to the free energy of a static quark anti-quark pair at infinite separation Foo{T) 



LreniT) = exp(-Foo(T)/(2r)), 
and obtained from the bare Polyakov loop as 

/l ^"-^ 

Lren{T) = z{Pf^Uare{P) = z{Pf^ ( -Tr J] f^o(^0,^) 

V xo=0 



(4) 



(5) 



Here the multiplicative renormalization constant z{j3) is related to the additive normalization 
of the potential c{j3) as z{l3) = exp(— c(/3)/2) discussed above. To make the comparison with 
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Figure 4. The renorinalized Polyakov loop 
as function of the temperature calculated for 
the HISQ, stout, asqtad and p4 actions at the 
physical value of the light quark mass. 



stout data the latter have to be multiplied by exp(— 0.91/(roT)), as in the stout calculation the 
potential was normalized to zero at r = tq. 

Our results for the quark mass mi = 0.05?7is on 32'^ x 8 lattices are shown in Fig. H] and 
compared to the stout, asqtad and p4 calculations. As one can see the HISQ calculations agree 
reasonably well with the stout results if the scale is set by rg in the stout calculations q On 
the other hand, the renormalized Polyakov loop calculated with the p4 and asqtad actions is 
noticeably smaller at low temperatures. At temperatures T > 200 MeV we see good agreement 
for different actions. The decrease of Foo(T), and thus the increase in the Polyakov loop could 
be related to onset of screening at high temperatures (e.g. see discussion in Ref. [27j). On 
the other hand, in the low-temperature region the increase of Lren is related to the fact that 
there are many static-light meson states that can contribute to the static quark free energy close 
to the transition temperature, while far away from the transition temperature it is determined 
by the binding energy of the lowest static-light mesons. The large taste symmetry breaking in 
the static-light meson sector is probably responsible for the discrepancy between p4 and HISQ 
results in the low-temperature region. Overall the differences between p4 and HISQ results are 
the smallest for the renormalized Polyakov loop as we will see in the following sections. This is 
partly due to the fact that the Polyakov loop is purely a gluonic observable. The taste symmetry 
breaking associated with light quarks enters only through loops, and the corresponding effect is 
smaller than for the quark number susceptibilities, for example. 

4. The chiral transition 

In the limit of zero light quark masses QCD has a chiral symmetry and the finite temperature 
transition is a true phase transition. The order parameter for this transition is the light chiral 
condensate {ip'>p)i- However, even at finite values of the quark mass the chiral condensate will 
show a rapid change in the transition region indicating an effective restoration of the chiral 
symmetry. Since the chiral condensate has an additive ultraviolet renormalization we consider 
the subtracted chiral condensate US] 



Ai,s(T) 






(6) 



Here the subscripts / and s refer to the light and strange quark condensates respectively 
normalized per single flavor, while the subscripts and r refer to the zero and finite temperature 



In what follows we will show the stout results using the temperature scale set by ro instead of the kaon decay 
constant /a-. We used the pubhshed values of ro and Jk in Refs. [181119] to convert the two scales. 
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Figure 5. The subtracted chiral condensate 
for the HISQ action compared with calcula- 
tions performed for the stout [19], asqtad [20] 
and p4 actions [17] . 



cases. In Fig. [5] the renormalized chiral condensate calculated with the HISQ action is shown 
and compared with results obtained with the stout [19], asqtad [20] and p4 [17] results. Our 
results agree reasonably well with the Nj- = 12 stout results. On the other hand, the subtracted 
chiral condensate is considerably smaller than for the asqtad and p4 actions on A'^^ = 8. This is 
due to the larger taste violating effects in the asqtad and p4 case. At finer lattice spacings the 
taste symmetry violations in the asqtad calculations become smaller and therefore the agreement 
between the asqtad calculations on A'^,- = 12 lattice and the HISQ calculations is much better 
(see Fig. [5]). We also note that if fx is used to set the scale there is no agreement between 
HISQ and stout calculations. 

To study the chiral aspects of the QCD transition we also consider the chiral susceptibility 
defined as the derivative of the light chiral condensate with respect to the quark masci 
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where Mi = D + 2mi is the staggered fermion matrix for light quarks. The first term in the above 
formula describes the fluctuation of the chiral condensate and is also called disconnected chiral 
susceptibility. The second term is called the connected chiral susceptibility and corresponds to 
the integrated scalar meson correlation function. For sufficiently small quark masses the chiral 
susceptibility is dominated by the disconnected part. The transition temperature is defined as 
the location of the peak of the chiral susceptibility. In Fig. [6] the disconnected and connected 
chiral susceptibilities for the HISQ action are displayed. We see that both of them show a peak- 
like structure for temperatures around 170 MeV. The disconnected chiral susceptibility has a 
broad shoulder below the transition temperature. This is due to the effect of the Goldstone 
modes [26]. In Fig. [7] we compare the HISQ results for the disconnected chiral susceptibility 
with the corresponding results for asqtad calculations performed at O.OSm^ [20]. There is a good 
agreement between A^,- = 8 HISQ and A',- = 12 asqtad calculations for the peak position. 

While the disconnected chiral susceptibility is finite in the continuum limit, the connected part 
is quadratically divergent for non-vanishing quark masses. To remove this additive ultraviolet 
divergences the Budapest- Wuppertal group subracted the zero temperature piece from the chiral 
susceptibility [18]. To remove the multiplicative renormalization they also multiplied the chiral 
condensate by the light quark mass squared [18], i.e. they considered the quantity 
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This is tiie 2-flavor chiral susceptibility as mi denotes the light quark mass which is common for u and d quarks. 
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Figure 6. The disconnected (left) and connected (right) chiral susceptibihties calculated for 
the HISQ action for the two values of the quark masses. 
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Figure 7. The comparison of the HISQ and 
asqtad [20j results for the disconnected chiral 
susceptibility at mi = 0.05ms. 



We calculated the same quantity for the HISQ action and show the comparison in Fig. [8l There 
is a reasonable agreement between HISQ and stout results if the tq is used to fix the scale. 
However, if the fx is used to set the scale the stout data shift to smaller temperatures. 



5. Fluctuations of conserved charges 

Fluctuations and correlations of conserved charges are good probes of deconfinement because 
they are sensitive to the underlying degrees of freedom, i.e. they can tell whether the relevant 
degrees of freedom of the system at a given temperature are hadronic or partonic. Here we 
consider quadratic fluctuations and correlations of conserved charges defined as 
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Here indices i and j refer to different conserved charges, like quark number, baryon number, 
strangeness, isospin etc. 
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Figure 8. The renormalized chiral susceptibility defined by Eq. ([8]) for the stout action with 
lattice spacing fixed by tq (left) and fx (right) [18]. We compare the stout results with the 
HISQ calculations at mi = 0.05?tt.s. 
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Figure 9. The strangeness fluctuations 
calculated with the HISQ action and compared 
to the asqtad action [20] for mi = 0.2m,s. 



We start our discussion with the strangeness fluctuation, since this quantity is often discussed 
in the literature in connection with deconfinement |161 [T8\ [19] . At low temperatures strangeness 
is carried by massive hadrons and therefore its fluctuations are suppressed. At high temperatures 
strangeness is carried by quarks and the effect of the non-zero strange quark mass is small. 
Therefore, in the transition region the strangeness fluctuation rises and eventually reaches 
a value close to that of an ideal quark gas. For mi = 0.2m,s our numerical results for the 
strangeness fluctuations are shown in Fig. [9] and compared with asqtad results. As one can see 
the strangeness fluctuations are larger for the HISQ action at low temperatures, T < 210 MeV. 
In other words, the transition region in the HISQ calculation shifts toward smaller temperatures. 
This behavior is in fact expected. Due to smaller taste symmetry violation pseudo-scalar masses 
as well as other hadron masses are smaller and therefore strangeness fluctuations are larger. For 
the physical quark mass, m/ = O.OSm^ we compare our calculations with the results obtained 
using the p4 action [T7] and the stout action [18l[19]. At low temperatures, T < 200 MeV the 
HISQ results are significantly larger than the p4 results but are in good agreement with the 
stout results. At high temperatures, T > 200 MeV the strangeness fluctuations calculated with 
the HISQ action are in reasonable agreement with the p4 results as well as the A'^,- = 12 stout 
results. We also see that in this temperature region the stout results show some cutoff (A,-) 
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Figure 10. The strangeness fluctuations calculated with the HISQ action for mi = 0.05m-s and 
compared to the results obtained with the p4, asqtad and stout actions. The lines show the 
prediction of the HRG model. 



dependence. This is due to the fact that the tree level ©(a^) lattice artifacts are not removed 
in the stout action. In the low temperature region we expect that the Hadron Resonance Gas 
(HRG) model gives a reasonably good description of the thermodynamic quantities, including 
strangeness fluctuations. Therefore in Fig. [10] we show the prediction of the HRG model. As one 
can see all lattice results fall below the HRG model result, although the difference between the 
lattice data and HRG is the smallest for the HISQ and stout actions. The discrepancy between 
the HRG model and the lattice data becomes larger at smaller temperatures, although the model 
is expected to be more reliable there. As explained in Ref. [28| this is due to taste symmetry 
violations leading to the distortion of the hadron spectrum, especially in the pseudo-scalar meson 
sector. Due to non-negligible splitting in the pseudo-scalar meson masses the contribution of 
kaons to strangeness fluctuations appears to be smaller than expected in the continuum [28]. 
As has been discussed in Ref. |28j this is mostly due to the discretization errors in the hadron 
spectrum. Taking into account the lattice spacing dependence of the hadron masses which enter 
the HRG calculations it is possible to get good agreement with lattice data l28j. Therefore in the 
figure we also show the HRG model calculations with hadron masses evaluated at lattice spacing 
corresponding to N^ = 8 and A^,- = 12 lattices using the formulas given in Ref. [28]. As one can 
see from the figure there is a very good agreement between this modified HRG calculation and 
lattice results. We note again, that for T < 160 MeV the splitting in the pseudo-scalar sector is 
the dominant source of the cutoff effects in the strangeness susceptibility. Finally we note that, 
when comparing the lattice results with the prediction of HRG we include all resonances up to 
2.5 GeV. Including resonances only up to 2.0 GeV will not change the strangeness fluctuations 
significantly, but further lowering the mass cutoff of the resonances will decrease the fluctuations 
for T > 160 MeV. 

We study the quark number and baryon number fluctuations which are sensitive to the light 
(non-strange) quark or hadron sector. The numerical results are shown in Fig. [TTJ Due to the 
much smaller quark mass the taste symmetry breaking effects are much more important in this 
sector. As the consequence the deviations from the HRG result are larger for these quantities 
and are significant even for HISQ calculations. However, the deviations between HRG and lattice 
results obtained with the HISQ action are much reduced compared to the previous calculations 
performed with the p4 action on N^ = 4 and 6 lattices at mi = O.Itti^ [30]. 

Next let us consider strangeness-baryon number and u- and s-quark number correlations. The 
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Figure 11. Fluctuation of the baryon number (left) and quark number (right) calculated with 
the HISQ action. Also shown as the solid line is the prediction of the HRG model. The lattice 
results are compared with the previous calculations performed with the p4 action [30]. 



numerical results for these quantities are shown in Fig. [121 At low temperatures strangeness- 
baryon number correlations are due to the presence of strange baryons. At high temperatures 
strange quarks are the relevant degrees of freedom and carry 1/3 unit of the baryon number and 
minus one unit of strangeness. As the result strangeness baryon number correlation approaches 
— 1/3 at high temperatures. The correlations of and u- and s-quark numbers at low temperatures 
are due to the presence of strange hadrons. Here both strange mesons and baryons contribute to 
this quantity. At high temperatures quarks are weakly interacting and the correlation between u- 
and s-quark numbers is very small. In fact, in perturbation theory it starts at order g^ log(fi') and 
numerically gives a result which is 10^ times smaller than the quark number susceptibility [29]. 
At low temperatures correlation of the baryon number and strangeness is very well described by 
HRG. The lattice results for u- and s-quark number correlations obtained with the HISQ action 
are not incompatible with HRG, but due to the large statistical errors it is difficult to get more 
precise conclusions. The correlations of conserved charges have also been calculated with the p4 
action, and the corresponding results are presented in Fig. [T2j As one can see from the figure 
the p4 results differ significantly from the HISQ results in the low temperature region. This is 
again due to large taste symmetry violations in the p4 calculations. 

Finally let us discuss the isospin fluctuations. The numerical results for the isospin 
fluctuations are shown in Fig. [13] and compared with the light quark number fluctuations. 
The latter are also compared with the strangeness fluctuations. The isospin fluctuations are 
very similar to the light quark number fluctuations at high temperatures. At low temperatures 
and in the transition region there are differences between the light quark number and isospin 
fluctuations. The strangeness and the light quark number fluctuations are also very different 
at all temperatures T < 260 MeV. Only above this temperature they become similar. One 
possible reason for the observed difference could be the different sensitivity of these quantities 
to the singular part of the free energy density. As discussed in Ref. [I6] the light quark number 
fluctuations may be sensitive to the singular part of the free energy density, while strangeness 
fluctuations are not sensitive to it. Isospin fluctuations are also not sensitive to the singular 
part of the free energy density as they are not coupled to the sigma fleld ^Jo It remains to 
be seen if the differences between the isospin and light quark number susceptibilities are due to 
the singular part of the free energy density. Current statistical errors do not allow to make firm 
conclusions. 



P.P. thanks F. Karsch and K. Redlich for the discussions on this point. 
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quark number (right) calculated with the HISQ action. Also shown is the prediction of the HRG 
model. The lattice results are compared with the previous calculations performed with the p4 
action [301. 
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Figure 13. The light quark number susceptibility compared to two times the isospin 
susceptibility (left) and to the strangeness susceptibility (right) for mi = 0.05ms- 



Correlations of conserved charges have also been calculated for tti/ = 0.2ms. The results are 
shown in Fig. [TH The results are similar to those obtained for the lighter quark mass. But due 
to the large quark mass the corresponding hadron masses are larger and the agreement with 
HRG is not that good. In the figure we also compare our numerical results with the results from 
previous calculations done with the p4 action [30]. There are significant differences between 
HISQ and p4 results also for this quark mass. Thus cutoff effects are the dominant source of the 
difference between p4 and HISQ results. In Fig. [T5]we show the light quark number fiuctuations 
in comparison with the isospin and strangeness fluctuations for m^i = 0.2m,s- As one can see from 
the figure the differences between the light quark number and strangeness fluctuations as well 
as light quark number and isospin fluctuations are smaller for this quark mass. This could be 
due to the fact that the differences in the masses of strange and non-strange hadrons are smaller 
for m,i = 0.2m,s, and also due to the reduced sensitivity of the light quark number fluctuations 
to the singular part. 

To summarize, in this section we have studied fluctuations and correlations of different 
conserved charges as function of temperature. We have found that at low temperatures, 
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Figure 14. Cross-correlation of the baryon number and strangeness (left) and and u- and s- 
quark number (right) calculated with the HISQ action for mi = 0.2ms and compared with p4 
results. Also shown is the prediction of the HRG model. 
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Figure 15. The light quark number susceptibility compared to two times the isospin 
susceptibility (left) and to the strangeness susceptibility (right) for m; = 0.2ms. 



T < 170 MeV, the fluctuations can be well understood in terms of the resonance gas model. 
To get a quantitative agreement with HRG it may be necessary to consider modifications of 
the hadron spectrum due to finite lattice spacing in the HRG calculations. This has been 
shown in some detail for the strangeness fluctuations. At high temperatures, T > 250 MeV, the 
fluctuations and correlations can be understood in terms of weakly interacting quark gas. In the 
transition region 170 MeV < T < 250 MeV there are no well defined quasi-particles which carry 
the conserved charges, i.e. it is difficult to define the relevant degrees of freedom. Thus, the 
deconfinement transition is a very gradual process if defined as the transition from hadronic to 
partonic degrees of freedom and it is impossible to associate a meaningful transition temperature 
with it. Some fluctuations may be sensitive to the singular part of the free energy density. In this 
case the inflection points will be related to the chiral transition temperature. In the literature 
(see e.g. Refs. |18lll9j ) attempts to define the deconfinement transition temperature as inflection 
points in the strangeness fluctuations and the renormalized Polyakov loop have been made. From 
the above discussion it is clear that the deconfinement transition temperature defined this way 
has little to do with the deconfinement temperature in the limit of large or infinite quark mass, 
which is related to the light modes of the Z{N) center symmetry breaking. Both strangeness 
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fluctuations and the renormalized Polyakov loop are dominated by the regular part of the free 
energy density. Thus, inflection points of these quantities need not to be related to the chiral 
transition temperature. In particular, these inflection points could be higher than the chiral 
transition temperature. This is shown in Fig. [TBI where we compare the disconnected chiral 
susceptibility with the derivatives of the renormalized Polyakov loop and Xs- 



6. Trace anomaly 

We also calculated the trace anomaly as function of temperature using the HISQ action for 
mi = 0.2ms and 0.05m,s. This quantity is often considered in lattice calculations as the pressure 
and other thermodynamic quantities are obtained by integrating over the trace anomaly to a 
given temperature. Previous calculations with the p4 and asqtad action gave a trace anomaly 
which was significantly below the HRG result. This deviation from the resonance gas complicated 
the use of lattice equation of state in hydrodynamic models since it was not clear how to 
implement consistent freezout at the hadronic stage of the evolution (see the discussion in Ref. 
[28]). In Fig. [T7] we show our results for the trace anomaly, e — 3p for the HISQ action and 
compare them with previous calculations performed with the asqtad and p4 actions. While 
at high temperatures (T > 250 MeV) all lattice data agree, there are noticeable differences in 
the low-temperature region. The lattice results obtained with the p4 and asqtad actions are 
lower than the HISQ results. This is due to the smaller breaking of the taste symmetry in the 
calculations performed with the HISQ action. We also compare the lattice calculations with the 
parametrization of the trace anomaly obtained from the matching the HRG at low temperature 
T < 180 MeV with the fit to the p4 lattice data at high temperatures for T > 250MeV. In the 
intermediate region the parametrization is constrained by the value of the entropy density at high 
temperatures. This parametrization gives an entropy density which is 5% below the ideal gas 
value at temperature of 800 MeV, and is called s95p-vl ^28J. It was also used in hydrodynamic 
models [28]. We also compare the lattice results with the parametrization obtained by Laine 
and Schroeder [H]. Within the statistical errors there is a good agreement between the lattice 
results obtained on A^,- = 8 lattice with m,i = 0.05ms and the s95p-vl parametrization. In the 
peak region the trace anomaly calculated on A^,- = 6 lattice is significantly below the A^,- = 8 
results, while at high temperatures the two calculations agree. Thus, the difference between 
A^T- = 6 and A^,- = 8 calculations follows the expected pattern of the cutoff dependence of the 
pressure in the free theory, namely at high temperatures the pressure obtained on A'^^ = 6 lattice 
should be below the pressure obtained on A',- = 8 lattice if the Naik term is used [32]. We also 
compare the lattice results with the HotQCD parametrization of the trace anomaly |16j . As can 
be seen in Fig. [17] the HotQCD parametrization agrees with s95p-vl parametrization at high 
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Figure 17. The trace anomaly calculated with the HISQ action for COSm^ (left) and 0.2ms 
(right) and compared with p4 and asqtad calculations. The open symbols in the left figure show 
the A'^^ = 6 results. We also compare the lattice results with the parametrizations of the trace 
anomaly: s95p-vl from Ref. |28] and the one by Laine and Schroeder |31j . In the right figure 
we also show the HotQCD parametrization of the trace anomaly |16j . 



temperatures. In the peak region and in the low-temperature region the HotQCD and s95p-vl 
parametrizations are different. This is due to the fact that the former is based on Nr = 8 p4 
and asqtad data in the entire temperature range. 



7. Conclusions 

In this paper we discussed recent results on QCD thermodynamics with the HISQ action. The 
use of this action significantly reduces the effect of the taste symmetry violation in the staggered 
fermion formulation. In particular, the quadratic splitting of pseudo-scalar mesons is the smallest 
for the HISQ action at a given value of the lattice spacing. The reduced lattice artifacts result 
in more realistic hadron spectrum. For instance, the vector meson masses are close to their 
continuum values for lattice spacings corresponding to the crossover region on A^,- = 8 lattices 
along the line of constant physics with (almost physical) mi = 0.05ms. We have calculated 
different thermodynamic quantities with the HISQ action on A'^^ = 6 and Nj- = 8 lattices, 
including the renormalized Polyakov loop, the chiral condensate, chiral susceptibility, trace 
anomaly and the fluctuations and correlations of conserved charges. We have found that at low 
temperatures the Polyakov loop agrees well with the results obtained with the stout action, if 
the lattice spacing is set by the Sommer scale ro, while at high temperatures it agrees with 
the p4 and asqtad calculations. The renormalized chiral condensate and chiral susceptibility 
obtained with the HISQ action agree quite well with the stout calculations but disagree with 
the results obtained with the p4 and asqtad actions on Nr = 8 lattices. However, the asqtad 
results on A^,- = 12 lattices agree well with the ones obtained with the HISQ action on A^,- = 8. 
We studied in detail different fluctuations and correlations of conserved charges. When 
comparing our results with previous ones obtained with the p4 and asqtad actions we have 
found differences in the low-temperature region. Those differences are due to the reduced taste 
symmetry violation in the HISQ action which gives a more realistic hadron spectrum. We have 
compared the lattice results in the low-temperature region with the prediction of the hadron 
resonance gas. We have found good agreement between the lattice results and the HRG once the 
cutoff effects have been taken into account. We have also found that the agreement between the 



lattice results and the physical HRG is best for the HISQ action. The deconfinenient transition 
can be understood as the transition from hadronic to quark degrees of freedom as seen by the 
different fluctuations/correlations. We have found that this transition is not abrupt but rather 
gradual. This implies that it is not possible to define a meaningful deconfinement transition 
temperature. 

Finally we considered the trace anomaly. The calculations performed with the HISQ action 
give results which agree quite well with the previous results obtained with the p4 and asqtad 
actions on A^,- = 6 and A^t- = 8 lattices at high temperatures. Differences are seen in the 
transition and in the low-temperature regions. The HISQ results on the trace anomaly agree 
well with the HRG in the low temperature region, and overall they agree with the s95p-vl 
parametrization that has been recently proposed for use in hydrodynamic models [28j. 
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